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We examine diffractive proton-proton scattering pp — > pX. Using a functional integral approach we derive 
the scattering amplitudes, which are governed by the expectation value of lightlike Wegner-Wilson loops. This 
expectation value is then evaluated using a cumulant expansion and the model of the stochastic vacuum. From 
the scattering amplitudes we calculate total and differential cross sections for high centre of mass energy and 
small momentum transfer and compare with experiments. 



1. Introduction 

In this article we study diffractive proton- 
proton scattering pp — > pX at high centre of mass 
(cm.) energies yfs > 20 GeV and small momen- 
tum transfer squared \t\ < 0.5 GeV 2 . The low 
momentum transfer implies that one has to ap- 
ply nonperturbative methods to investigate these 
processes. 

A description of soft hadronic high energy re- 
actions, starting from a microscopic level, was 
developed in [Q where in the case of an abelian 
gluon model the pomeron properties were related 
to nonperturbative aspects of the vacuum like the 
gluon condensate introduced by Shifman, Vain- 
shtein and Zakharov J2|. These methods were 
generalized to QCD in ||]. The quantity gov- 
erning the scattering amplitude was found to be 
a correlation function of two lightlike Wegner- 
Wilson loops [^J. These correlation functions are 
evaluated using the model of the stochastic vac- 
uum (MSV) Q after an analytic continuation 
from Euclidian space to Minkowski space Q . 

In this paper we will apply the model to in- 
elastic diffractive proton-proton scattering. The 
hadronic scattering amplitude as derived in pLE, 
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0J|] is given in section |[ the result for the loop- 
loop correlation function is shown in section 0. 
The numerical results are presented in sectionH 
and compared to experimental data. 

2. The scattering amplitudes 

In this section we present our basic formulae for 
hadron-hadron scattering. Consider the reaction 



/ii(Pi) + h 2 {P 2 ) -» /ii(P 3 ) + X(P 4 ), 



(1) 



where h\ and h 2 are hadrons, X is again h 2 or a 
diffractive excitation of h 2 . The hadrons hi, h 2 
are modeled as qq and quark-diquark wave pack- 
ets for mesons and baryons respectively. For the 
wave functions we have chosen a Bauer-Stech- 
Wirbel type ansatz || . The diffractive final state 
X is modeled by a qq-pair (or quark-diquark pair) 
in a colour singlet state. For the description of 
the quark and antiquark (diquark), we use free 
plane waves. Integration over all allowed values 
in phase space and the closure relation then yield 
all possible diffractive final states X, where the 
case of elastic scattering is also included. 

In the framework of the model presented in 
we obtain the scattering amplitude for reac- 
tion (m) as 



Tfi = (2is)(2n) / db T bTJo(V^ibT)JdiS. (2) 
Jo 

Here Jo is the Bessel function of zeroth degree and 
Jdiff is the diffractive profile function for which we 
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obtain 

Jdiff(b T ,^') = - / d 2 x T d 2 y T J dzw 3 i(x T ,z) 

V2^V22 / (l-z')e _iA4I " yr ^a(yr,*0 
(w+(^b T + (i-z)x T ,x T ) 

W_(-^b T + (i-^)yr,yr)-l) G , (3) 

where 1031 (xy, z) denotes the profile function for 
the overlap between initial and final state of the 
hadron hi for fixed xt and z. ^(yT,^') defines 
the initial state wave function of h 2 . Using the 
ansatz of 11 we have 

v/2z(l-z)e-^-^ 2 / 4 <, (4) 

Wij(XT,z) = ipi(x. T ,z)cpj(x.T,z)- (5) 

Here 1^ is a normalization factor. The lightlikc 
Wegner- Wilson loops W± are given by 

W± :=±trPexp(-i fl ^ d^G£(z)y), (6) 

where P denotes path ordering and C± is the 
curve consisting of two lightlike worldlines for the 
quark and the antiquark (diquark) and connect- 
ing pieces at ±00. xj- and y-r define the exten- 
sion and orientation in transverse position space 
of the two loops representing the two hadrons hi 
and /12 respectively, z (z') parametrizes the frac- 
tion of the longitudinal momentum of hadron hi 
(/12) carried by the quark. The impact parameter 
is given by by. 

The symbol (. . .)g denotes the functional inte- 
gration which correlates the two loops. In (|^) the 
loop-loop correlation function is multiplied with 
the profile function W31, the incoming wave func- 
tion if2 and then integrated over all extensions 
and orientations of the loops in transverse space 
as well as over the longitudinal momentum frac- 
tion z of hadron hi , which is not diffractively ex- 
cited. The hadronic scattering amplitude is then 
given by a Fourier-Bessel transform of this expres- 
sion with respect to the impact parameter b^. 



3. Evaluation of the scattering amplitudes 

Now we perform the functional integral in (^) 
making use of the MSV. A detailed presenta- 
tion of the MSV can be found in [|]-||, where 
both the original formulation in Euclidian space- 
time and the analytic continuation to Minkowski 
space-time are discussed. 

The first step is to make a cumulant expansion |tJ 
for the loop-loop correlation function. Proceed- 
ing as explained in [ flO| and assuming \\\ <C 1 we 
find 

{W + W--1) G = -\ X \ (7) 



which is the result of the traditional expansion 
method 01 . Here we have defined 




where G2 is the gluon condensate, A = 8a/37r 
is connected to the correlation length a and k is 
a parameter which is related to the non-abclian 
character of the correlator. -^2,3 are the modified 
Bessel functions of second and third degree. The 
vectors r,-j with i = x,y and j = q,q run from the 
coordinate origin to the positions of the quarks 
and antiquarks in transverse space. 

4. Total and differential cross sections 

Now we consider diffractive proton-proton scat- 
tering, i.e. hi = h 2 = p. To calculate cross sec- 
tions, we have to fix our free parameters, those 
of the MSV and the parameters of the proton 
wave function, i.e. the extension parameter S p 
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Figure 1. The integrated single diffractive cross 
section as a function of yfs 
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Figure 2. The differential diffractive cross section 
da sd /dt [mb/GeV 2 ] at y/s = 23.5 GeV 



and the width of the longitudinal momentum dis- 
tribution z p . The set of MSV parameters used 
in this work has been established in JD]]: G2 = 
(501 MeV) 4 , a = 0.346 fm, n = 0.74. The pro- 
ton extension parameter is allowed to be energy 
dependent. From a fit of the total cross section 
as calculated in the model to the pomcron part 
of the Donnachic-Landshoff (DL) parametrisation 
for otot [f2| we obtain the following connection 



between S p and s: 



S p (s) = 0.624 



(— 

VGeV 2 



fm. 



(10) 



The width of the longitudinal momentum distri- 
bution has been chosen as z p = 0.21 which gives 
a best fit to the electric form factor of the proton 
calculated in the framework of our model fli3fl . 

With all parameters fixed, we can now do the 
numerical calculations for the scattering ampli- 
tude. We calculate the differential diffractive 
cross section from ^ , where we insert (^) for the 
loop-loop correlation function and the wave and 
profile functions (|j), (||) in (|^) 



do-diff 



where 



d b V = 



(2nf±-\T fl \ 2 d 5 V, 



(2tt) 9 4s z'{1 - z 1 ) 



d 2 P 4T d 2 A AT dz' 



(11) 



(12) 



is the 5-dimensional phase space measure. 

To obtain the cross section for single diffrac- 
tive dissociation, we have to subtract the elastic 
contribution (cf. |10[] ) and then multiply by 2 to 
account for the reaction where the first proton 
breaks up. We then find for the integrated single 
diffractive cross section as a function of y/s the 
result shown in Fig. [l], where we have included 
our result for the integrated elastic cross section 
for comparison. Comparing our results to exper- 
imental data, one has to keep in mind that the 
overall normalisation uncertainty of the experi- 
ments is of O(10%). 

We now calculate the fraction R — a s< i/((j e \ + 
cr s d) of the diffractive dissociation cross section 
to the sum of the diffractive dissociation and the 
elastic parts and compare to experiment 1 19 l4| - 
jL8| . For cm. energies y/s = 23.5, 546 and 
1800 GeV we find R = 0.49, 0.46 and 0.45, re- 
spectively, compared to the experimental values 
R = 0.49±0.01, 0.41±0.02 and 0.38±0.02. As can 
be seen from Fig. [I] and the calculated i?-values, 
our model is in qualitative agreement with the ex- 
perimentally observed trend that the diffractive 
dissociation cross section grows more slowly with 
increasing energy than the elastic cross section. 

The result for the differential cross section of 
the diffractive dissociation is shown in Fig. ||. A 
fit to our result of the form do~ s d/dt = A expat 
gives b — 7.3±0.2 GeV -2 compared to the exper- 
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imental value b = 7.0 ± 0.3 GeV~ 2 @, meaning 
that our model can reproduce the shape of the 
diffractive |t|-distribution reasonably well in the 
range of small |i| considered here. 

5. Conclusion 

In this work we have calculated integrated and 
differential cross sections for inelastic diffractive 
scattering at high cm energies and small momen- 
tum transfer from the scattering amplitudes. In 
the ISR energy regime, our model describes the 
experimental data on proton-proton scattering 
well within the numerical and experimental er- 
rors and reproduces the experimentally observed 
features of diffractive scattering (e.g. the shrink- 
age of the ratio of the cross sections of diffrac- 
tively excited to the sum of elastic scattering and 
diffractive excitation). Future investigations will 
cover the shapes of the diffractive mass spectra 
and the influence of higher terms in the cumulant 
expansion of the loop-loop correlation function. 

A further field of application for our model will 
be the study of double diffractive dissociation re- 
actions, where C = — 1 contributions (odderon 
exchange) would be included already in our cur- 
rent approximation. Upcoming experiments at 
RHIC will be a rich source for new experimen- 
tal data for both single and double diffractive re- 
actions in hadronic reactions at high cm. ener- 
gies and therefore the study of inelastic diffractive 
scattering will remain an interesting and instruc- 
tive subject for future studies. 
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